This analysis document compliments FIA NLS Models: Biomass vs. Stand Age. All of the background information from that document applies to these analyses, which are extensions to them. The difference between that document and this analysis is the use of different data subsets.
Here, we fit the models using: 1) a temporally-balanced dataset, where we take the first and most-recent plot record for all plots in the dataset, 2) a temporally-balanced dataset (same as #1), but which excludes plot locations which have experienced harvest (at any point over the study interval 2000-2022)
Below the model fitting procedure is implemented by ecoprovince:
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 4834 1203.6
## 2 4833 1092.6 1 111.05 491.22 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 52557.62
## 2 2 52091.41
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.45951 0.16992 2.704 0.00687 **
## alpha 0.84855 0.03503 24.221 < 2e-16 ***
## A 473.76484 37.00356 12.803 < 2e-16 ***
## k 205.69316 18.21698 11.291 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4755 on 4833 degrees of freedom
##
## Number of iterations to convergence: 7
## Achieved convergence tolerance: 2.017e-06
## (1 observation deleted due to missingness)
## Error in Mod.Sel3 %in% c(1, "1a", "1b", "1c", 4) :
## object 'Mod.Sel3' not found
## Error in Mod.Sel3 %in% c(1, "1a", "1b", "1c", 4) :
## object 'Mod.Sel3' not found
## Error in Mod.Sel3 %in% c(1, "1a", "1b", "1c", 4) :
## object 'Mod.Sel3' not found
## model AIC
## 1 2 52091.41
## 2 2a NA
## 3 2b NA
## 4 2c NA
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.45951 0.16992 2.704 0.00687 **
## alpha 0.84855 0.03503 24.221 < 2e-16 ***
## A 473.76484 37.00356 12.803 < 2e-16 ***
## k 205.69316 18.21698 11.291 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4755 on 4833 degrees of freedom
##
## Number of iterations to convergence: 7
## Achieved convergence tolerance: 2.017e-06
## (1 observation deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 4832 1179.6
## 2 4831 1067.7 1 111.87 506.17 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 52091.41
## 2 4 52463.90
## 3 5 51983.93
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.58118 0.17752 3.274 0.00107 **
## alpha 0.84078 0.03411 24.646 < 2e-16 ***
## a 33.40070 1.99286 16.760 < 2e-16 ***
## b 118.88478 5.92661 20.059 < 2e-16 ***
## c 118.37512 5.85956 20.202 < 2e-16 ***
## d 1.01653 0.05233 19.425 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4701 on 4831 degrees of freedom
##
## Algorithm "port", convergence message: relative convergence (4)
## (1 observation deleted due to missingness)
## Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
## ℹ Please use `linewidth` instead.
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 12943 5165.4
## 2 12942 5044.1 1 121.33 311.3 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 135621.4
## 2 2 135315.7
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.95237 0.12759 7.464 8.91e-14 ***
## alpha 0.62169 0.03318 18.736 < 2e-16 ***
## A 165.96924 5.63575 29.449 < 2e-16 ***
## k 78.87707 3.15407 25.008 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.6243 on 12942 degrees of freedom
##
## Number of iterations to convergence: 4
## Achieved convergence tolerance: 8.623e-06
## (16 observations deleted due to missingness)
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 12942 5044.1
## 2 12941 5021.8 1 22.266 57.378 3.839e-14 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 135315.7
## 2 2a 135260.4
## 3 2b 135317.6
## 4 2c NA
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 9.238e-01 1.259e-01 7.336 2.33e-13 ***
## alpha 6.214e-01 3.094e-02 20.082 < 2e-16 ***
## A 1.921e+02 9.346e+00 20.559 < 2e-16 ***
## k 1.085e+02 7.979e+00 13.593 < 2e-16 ***
## p 2.121e-02 2.492e-03 8.509 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.6229 on 12941 degrees of freedom
##
## Number of iterations to convergence: 9
## Achieved convergence tolerance: 5.063e-06
## (16 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 12941 5088.9
## 2 12940 4912.4 1 176.51 464.95 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2a 135260.4
## 2 4 135432.3
## 3 5 134977.3
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.19051 0.13607 8.749 <2e-16 ***
## alpha 0.67333 0.02702 24.922 <2e-16 ***
## a 13.38590 0.51341 26.073 <2e-16 ***
## b 75.80154 2.36631 32.034 <2e-16 ***
## c 121.05709 5.45544 22.190 <2e-16 ***
## d 1.41595 0.04287 33.030 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.6161 on 12940 degrees of freedom
##
## Algorithm "port", convergence message: relative convergence (4)
## (16 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 5442 948.80
## 2 5441 857.87 1 90.926 576.69 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 59550.63
## 2 2 59004.09
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.1644 0.1261 1.304 0.192
## alpha 0.8050 0.0312 25.799 <2e-16 ***
## A 500.5614 30.8115 16.246 <2e-16 ***
## k 155.6108 11.1725 13.928 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.3971 on 5441 degrees of freedom
##
## Number of iterations to convergence: 7
## Achieved convergence tolerance: 2.272e-06
## (1 observation deleted due to missingness)
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 5441 857.87
## 2 5440 856.53 1 1.3429 8.529 0.00351 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 59004.09
## 2 2a 58997.56
## 3 2b 58973.93
## 4 2c NA
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.23383 0.13009 1.797 0.0723 .
## alpha 0.81624 0.03142 25.975 <2e-16 ***
## A 293.53650 21.44552 13.688 <2e-16 ***
## k 64.20469 6.41740 10.005 <2e-16 ***
## s 1.31952 0.05962 22.132 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.3959 on 5440 degrees of freedom
##
## Number of iterations to convergence: 7
## Achieved convergence tolerance: 2.275e-06
## (1 observation deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 5440 940.37
## 2 5439 848.31 1 92.062 590.26 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2b 58973.93
## 2 4 59506.03
## 3 5 58947.03
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.24222 0.13018 1.861 0.0629 .
## alpha 0.80822 0.03087 26.178 < 2e-16 ***
## a 21.08256 2.67854 7.871 4.22e-15 ***
## b 179.51639 10.45102 17.177 < 2e-16 ***
## c 156.83930 15.71326 9.981 < 2e-16 ***
## d 1.51879 0.08979 16.915 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.3949 on 5439 degrees of freedom
##
## Algorithm "port", convergence message: relative convergence (4)
## (1 observation deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 3547 1178.6
## 2 3546 1136.5 1 42.068 131.25 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 38157.03
## 2 2 38030.00
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.03509 0.24350 4.251 2.18e-05 ***
## alpha 0.72686 0.05912 12.295 < 2e-16 ***
## A 445.56009 55.43529 8.037 1.24e-15 ***
## k 254.71684 34.82136 7.315 3.17e-13 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5661 on 3546 degrees of freedom
##
## Number of iterations to convergence: 8
## Achieved convergence tolerance: 2.413e-06
## (2 observations deleted due to missingness)
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 3546 1136.5
## 2 3545 1132.9 1 3.5914 11.238 0.00081 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 38030.00
## 2 2a 38020.77
## 3 2b NA
## 4 2c NA
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.091545 0.248065 4.400 1.11e-05 ***
## alpha 0.753308 0.059393 12.683 < 2e-16 ***
## A 342.271688 43.589966 7.852 5.38e-15 ***
## k 170.045425 28.883552 5.887 4.29e-09 ***
## p -0.015140 0.006514 -2.324 0.0202 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5653 on 3545 degrees of freedom
##
## Number of iterations to convergence: 10
## Achieved convergence tolerance: 5.957e-06
## (2 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 3545 1160.1
## 2 3544 1108.8 1 51.282 163.91 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2a 38020.77
## 2 4 38104.80
## 3 5 37946.29
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.21145 0.25520 4.747 2.14e-06 ***
## alpha 0.76881 0.05439 14.134 < 2e-16 ***
## a 14.01898 1.39169 10.073 < 2e-16 ***
## b 101.92211 6.02465 16.918 < 2e-16 ***
## c 115.56990 8.21593 14.067 < 2e-16 ***
## d 1.19701 0.06667 17.955 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5593 on 3544 degrees of freedom
##
## Algorithm "port", convergence message: relative convergence (4)
## (2 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 6383 1208.6
## 2 6382 1143.1 1 65.521 365.81 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 66739.63
## 2 2 66385.70
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.80821 0.12030 6.718 2e-11 ***
## alpha 0.70093 0.03432 20.426 <2e-16 ***
## A 212.97304 8.46135 25.170 <2e-16 ***
## k 72.75583 3.94421 18.446 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4232 on 6382 degrees of freedom
##
## Number of iterations to convergence: 5
## Achieved convergence tolerance: 2.25e-06
## (2 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 6382 1143.1
## 2 6381 1142.5 1 0.6350 3.5469 0.0597 .
## 3 6381 1136.2 0 0.0000
## 4 6380 1128.1 1 8.1539 46.1159 1.215e-11 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 66385.70
## 2 2a 66384.15
## 3 2b 66349.16
## 4 2c 66305.16
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.85759 0.12173 7.045 2.05e-12 ***
## alpha 0.71209 0.03381 21.062 < 2e-16 ***
## A 127.11413 4.61707 27.531 < 2e-16 ***
## k 37.02397 1.25609 29.475 < 2e-16 ***
## p 0.14094 0.01883 7.486 8.09e-14 ***
## s 2.17271 0.15279 14.220 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4205 on 6380 degrees of freedom
##
## Number of iterations to convergence: 5
## Achieved convergence tolerance: 4.063e-06
## (2 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 6381 1195.9
## 2 6380 1127.7 1 68.139 385.5 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2c 66305.16
## 2 4 66675.81
## 3 5 66303.16
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.84613 0.12109 6.988 3.08e-12 ***
## alpha 0.71122 0.03383 21.023 < 2e-16 ***
## a 17.80094 2.07031 8.598 < 2e-16 ***
## b 98.51432 3.90582 25.222 < 2e-16 ***
## c 113.18502 6.98395 16.206 < 2e-16 ***
## d 1.44570 0.07453 19.398 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4204 on 6380 degrees of freedom
##
## Algorithm "port", convergence message: relative convergence (4)
## (2 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 7773 2698.5
## 2 7772 2568.4 1 130.03 393.48 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 85845.56
## 2 2 85463.52
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.24476 0.17019 7.314 2.86e-13 ***
## alpha 0.57095 0.02662 21.444 < 2e-16 ***
## A 257.12975 11.27279 22.810 < 2e-16 ***
## k 67.16853 3.01370 22.288 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5749 on 7772 degrees of freedom
##
## Number of iterations to convergence: 6
## Achieved convergence tolerance: 1.572e-06
## (14 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 7772 2568.4
## 2 7771 2545.5 1 22.931 70.004 < 2.2e-16 ***
## 3 7771 2566.4 0 0.000
## 4 7770 2496.6 1 69.841 217.364 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 85463.52
## 2 2a 85395.79
## 3 2b 85459.40
## 4 2c 85246.85
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.506e+00 1.836e-01 8.201 2.76e-16 ***
## alpha 6.512e-01 2.258e-02 28.840 < 2e-16 ***
## A 1.450e+02 6.248e+00 23.205 < 2e-16 ***
## k 3.016e+01 1.222e+00 24.686 < 2e-16 ***
## p 1.076e-01 7.674e-03 14.024 < 2e-16 ***
## s 1.963e+00 9.744e-02 20.143 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5668 on 7770 degrees of freedom
##
## Number of iterations to convergence: 12
## Achieved convergence tolerance: 5.033e-06
## (14 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 7771 2669.8
## 2 7770 2494.5 1 175.36 546.22 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2c 85246.85
## 2 4 85766.57
## 3 5 85240.29
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.47599 0.18195 8.112 5.75e-16 ***
## alpha 0.65085 0.02255 28.860 < 2e-16 ***
## a 16.12785 0.82771 19.485 < 2e-16 ***
## b 120.36292 5.58045 21.569 < 2e-16 ***
## c 112.83900 9.26790 12.175 < 2e-16 ***
## d 1.63350 0.06625 24.657 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5666 on 7770 degrees of freedom
##
## Algorithm "port", convergence message: relative convergence (4)
## (14 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 7905 4072.2
## 2 7904 3928.6 1 143.61 288.94 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 88878.43
## 2 2 88596.51
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.76811 0.18004 4.266 2.01e-05 ***
## alpha 0.55314 0.02987 18.521 < 2e-16 ***
## A 270.07655 14.26944 18.927 < 2e-16 ***
## k 71.69928 3.93397 18.226 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.705 on 7904 degrees of freedom
##
## Number of iterations to convergence: 6
## Achieved convergence tolerance: 1.137e-06
## (32 observations deleted due to missingness)
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 7904 3928.6
## 2 7903 3859.0 1 69.597 142.53 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 88596.51
## 2 2a 88457.16
## 3 2b 88567.74
## 4 2c NA
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 9.272e-01 1.885e-01 4.920 8.85e-07 ***
## alpha 6.432e-01 2.456e-02 26.191 < 2e-16 ***
## A 4.083e+02 4.091e+01 9.981 < 2e-16 ***
## k 1.596e+02 2.131e+01 7.489 7.67e-14 ***
## p 2.506e-02 1.901e-03 13.185 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.6988 on 7903 degrees of freedom
##
## Number of iterations to convergence: 7
## Achieved convergence tolerance: 4.089e-06
## (32 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 7903 4029.5
## 2 7902 3773.8 1 255.69 535.39 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2a 88457.16
## 2 4 88799.14
## 3 5 88282.71
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.00852 0.19081 5.285 1.29e-07 ***
## alpha 0.71043 0.02138 33.233 < 2e-16 ***
## a 20.51866 1.02371 20.043 < 2e-16 ***
## b 121.27486 6.30094 19.247 < 2e-16 ***
## c 111.11452 9.25971 12.000 < 2e-16 ***
## d 1.50283 0.06832 21.998 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.6911 on 7902 degrees of freedom
##
## Algorithm "port", convergence message: relative convergence (4)
## (32 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 826 218.04
## 2 825 194.59 1 23.448 99.412 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 9094.121
## 2 2 9001.802
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 3.865e-01 4.112e-01 0.940 0.3476
## alpha 7.628e-01 6.984e-02 10.922 <2e-16 ***
## A 1.255e+03 5.318e+02 2.360 0.0185 *
## k 4.893e+02 2.232e+02 2.192 0.0286 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4857 on 825 degrees of freedom
##
## Number of iterations to convergence: 7
## Achieved convergence tolerance: 1.631e-06
## (1 observation deleted due to missingness)
## Error in nls(get(paste("f_", Mod.Sel1, "a", sep = "")), data = G_234, :
## step factor 0.000488281 reduced below 'minFactor' of 0.000976562
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 825 194.59
## 2 823 192.48 2 2.1059 4.5022 0.01136 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 9001.802
## 2 2a NA
## 3 2b NA
## 4 2c 8996.781
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.39536 0.41113 0.962 0.337
## alpha 0.78925 0.07021 11.242 < 2e-16 ***
## A 395.47217 189.29921 2.089 0.037 *
## k 99.96670 70.69953 1.414 0.158
## p -0.01220 0.01423 -0.857 0.392
## s 1.13109 0.24438 4.628 4.28e-06 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4836 on 823 degrees of freedom
##
## Number of iterations to convergence: 40
## Achieved convergence tolerance: 3.1e-06
## (1 observation deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 824 216.69
## 2 823 192.29 1 24.399 104.43 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2c 8996.781
## 2 4 9092.966
## 3 5 8995.935
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.4024 0.4120 0.977 0.32908
## alpha 0.7883 0.0701 11.246 < 2e-16 ***
## a 0.0000 5.9324 0.000 1.00000
## b 340.8238 204.3728 1.668 0.09576 .
## c 675.4071 940.5067 0.718 0.47288
## d 2.5547 0.6884 3.711 0.00022 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4834 on 823 degrees of freedom
##
## Algorithm "port", convergence message: relative convergence (4)
## (1 observation deleted due to missingness)
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 1387 347.56
## 2 1386 336.63 1 10.934 45.02 2.834e-11 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 14526.30
## 2 2 14483.87
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.28977 0.25193 1.150 0.25
## alpha 0.68210 0.09544 7.147 1.43e-12 ***
## A 261.45525 28.66195 9.122 < 2e-16 ***
## k 106.87482 14.79744 7.223 8.39e-13 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4928 on 1386 degrees of freedom
##
## Number of iterations to convergence: 5
## Achieved convergence tolerance: 9.101e-06
## (2 observations deleted due to missingness)
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 1386 336.63
## 2 1385 336.63 1 5.7789e-05 2e-04 0.9877
## model AIC
## 1 2 14483.87
## 2 2a 14485.87
## 3 2b 14469.11
## 4 2c NA
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.34766 0.25540 1.361 0.174
## alpha 0.70227 0.09501 7.391 2.51e-13 ***
## A 142.90894 13.48623 10.597 < 2e-16 ***
## k 39.02880 3.98802 9.787 < 2e-16 ***
## s 1.74387 0.18179 9.593 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.49 on 1385 degrees of freedom
##
## Number of iterations to convergence: 8
## Achieved convergence tolerance: 6.83e-06
## (2 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 1385 337.23
## 2 1384 325.82 1 11.409 48.464 5.181e-12 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2b 14469.11
## 2 4 14488.37
## 3 5 14442.53
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.34162 0.25201 1.356 0.175
## alpha 0.70124 0.09431 7.435 1.82e-13 ***
## a 24.39435 4.05270 6.019 2.24e-09 ***
## b 93.28179 7.41989 12.572 < 2e-16 ***
## c 102.73710 8.92197 11.515 < 2e-16 ***
## d 1.13595 0.10865 10.455 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4852 on 1384 degrees of freedom
##
## Algorithm "port", convergence message: relative convergence (4)
## (2 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 439 237.09
## 2 438 235.09 1 2.0089 3.7429 0.05368 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 4609.958
## 2 2 4608.197
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau -0.1340 0.6019 -0.223 0.82393
## alpha 0.3504 0.1744 2.009 0.04512 *
## A 320.3752 106.7015 3.003 0.00283 **
## k 143.2259 52.7209 2.717 0.00686 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.7326 on 438 degrees of freedom
##
## Number of iterations to convergence: 8
## Achieved convergence tolerance: 5.634e-06
## (2 observations deleted due to missingness)
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 438 235.09
## 2 437 234.74 1 0.3499 0.6514 0.4201
## model AIC
## 1 2 4608.197
## 2 2a 4609.539
## 3 2b NA
## 4 2c NA
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau -0.1340 0.6019 -0.223 0.82393
## alpha 0.3504 0.1744 2.009 0.04512 *
## A 320.3752 106.7015 3.003 0.00283 **
## k 143.2259 52.7209 2.717 0.00686 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.7326 on 438 degrees of freedom
##
## Number of iterations to convergence: 8
## Achieved convergence tolerance: 5.634e-06
## (2 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 437 226.64
## 2 436 221.69 1 4.9467 9.7286 0.001934 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 4608.197
## 2 4 4594.028
## 3 5 4586.274
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau -0.2440 0.5501 -0.443 0.657650
## alpha 0.5283 0.1577 3.350 0.000877 ***
## a 9.5645 1.9483 4.909 1.30e-06 ***
## b 88.5138 12.4743 7.096 5.24e-12 ***
## c 55.9913 5.9220 9.455 < 2e-16 ***
## d 1.0392 0.1138 9.129 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.7131 on 436 degrees of freedom
##
## Algorithm "port", convergence message: relative convergence (4)
## (2 observations deleted due to missingness)
add p model: does not fit
add s model: does not fit
add s+p model: does not fit
unable to fit model (only 64 observations)
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
add p model: does not fit
add s model: does not fit
add s+p model: does not fit
unable to fit model (0 observations)
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## Error in if (any(nEQ <- vNms != make.names(vNms))) vNms[nEQ] <- paste0("`", :
## missing value where TRUE/FALSE needed
## Error in if (any(nEQ <- vNms != make.names(vNms))) vNms[nEQ] <- paste0("`", :
## missing value where TRUE/FALSE needed
## model AIC
## 1 1 NA
## 2 2 NA
## Warning in min(AIC1_322$AIC, na.rm = T): no non-missing arguments to min;
## returning Inf
## Error in get(paste("nls_322.", Mod.Sel1, sep = "")) :
## object 'nls_322.' not found
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 149 87.775
## 2 148 83.201 1 4.5741 8.1365 0.00496 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 1640.707
## 2 2 1634.572
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 3.0055 4.1286 0.728 0.46779
## alpha 0.9714 0.2977 3.263 0.00137 **
## A 231.9254 226.9888 1.022 0.30857
## k 219.8242 199.9830 1.099 0.27346
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.7498 on 148 degrees of freedom
##
## Number of iterations to convergence: 6
## Achieved convergence tolerance: 4.357e-06
## (2 observations deleted due to missingness)
## Error in nls(get(paste("f_", Mod.Sel1, "a", sep = "")), data = G_332, :
## number of iterations exceeded maximum of 50
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 148 83.201
## 2 147 83.108 1 0.093179 0.1648 0.6854
## model AIC
## 1 2 1634.572
## 2 2a NA
## 3 2b 1636.402
## 4 2c NA
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 3.0055 4.1286 0.728 0.46779
## alpha 0.9714 0.2977 3.263 0.00137 **
## A 231.9254 226.9888 1.022 0.30857
## k 219.8242 199.9830 1.099 0.27346
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.7498 on 148 degrees of freedom
##
## Number of iterations to convergence: 6
## Achieved convergence tolerance: 4.357e-06
## (2 observations deleted due to missingness)
## Error in nls(f_5, data = G_332, start = c(tau = tau.start, alpha = alpha.start, :
## Convergence failure: iteration limit reached without convergence (10)
## model AIC
## 1 2 1634.572
## 2 4 1636.879
## 3 5 NA
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 3.0055 4.1286 0.728 0.46779
## alpha 0.9714 0.2977 3.263 0.00137 **
## A 231.9254 226.9888 1.022 0.30857
## k 219.8242 199.9830 1.099 0.27346
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.7498 on 148 degrees of freedom
##
## Number of iterations to convergence: 6
## Achieved convergence tolerance: 4.357e-06
## (2 observations deleted due to missingness)
simple log-normal model: does not fit
log-normal phi model: does not fit
model not fitted because only 62 observations
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 5103 963.25
## 2 5102 841.07 1 122.19 741.21 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 54262.24
## 2 2 53571.63
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.63454 0.16362 3.878 0.000107 ***
## alpha 0.80068 0.02686 29.811 < 2e-16 ***
## A 446.29921 29.09383 15.340 < 2e-16 ***
## k 198.11601 13.75314 14.405 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.406 on 5102 degrees of freedom
##
## Number of iterations to convergence: 7
## Achieved convergence tolerance: 3.889e-06
## (2 observations deleted due to missingness)
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 5102 841.07
## 2 5101 833.07 1 8.0008 48.99 2.905e-12 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 53571.63
## 2 2a 53524.82
## 3 2b NA
## 4 2c NA
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.51537 0.15576 3.309 0.000944 ***
## alpha 0.81328 0.02703 30.083 < 2e-16 ***
## A 324.77328 20.26162 16.029 < 2e-16 ***
## k 102.66625 10.74934 9.551 < 2e-16 ***
## p -0.05130 0.01127 -4.551 5.47e-06 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4041 on 5101 degrees of freedom
##
## Number of iterations to convergence: 6
## Achieved convergence tolerance: 2.641e-06
## (2 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 5101 952.58
## 2 5100 828.03 1 124.55 767.12 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2a 53524.82
## 2 4 54209.36
## 3 5 53495.89
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.53028 0.15658 3.387 0.000713 ***
## alpha 0.81390 0.02671 30.470 < 2e-16 ***
## a 14.56503 2.98921 4.873 1.14e-06 ***
## b 155.65092 10.10200 15.408 < 2e-16 ***
## c 179.79825 19.79447 9.083 < 2e-16 ***
## d 1.58647 0.10396 15.261 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4029 on 5100 degrees of freedom
##
## Algorithm "port", convergence message: relative convergence (4)
## (2 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 5180 879.40
## 2 5179 824.61 1 54.79 344.11 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 57351.18
## 2 2 57019.77
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.69968 0.13285 5.267 1.45e-07 ***
## alpha 0.82371 0.04201 19.609 < 2e-16 ***
## A 266.46784 10.04611 26.524 < 2e-16 ***
## k 61.97132 3.14535 19.703 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.399 on 5179 degrees of freedom
##
## Number of iterations to convergence: 7
## Achieved convergence tolerance: 3.081e-06
## (3 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 5179 824.61
## 2 5178 823.80 1 0.8171 5.1361 0.02347 *
## 3 5178 817.30 0 0.0000
## 4 5177 805.02 1 12.2794 78.9674 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 57019.77
## 2 2a 57016.63
## 3 2b 56975.59
## 4 2c 56899.12
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.87739 0.14139 6.206 5.87e-10 ***
## alpha 0.82469 0.03959 20.832 < 2e-16 ***
## A 162.45200 5.61416 28.936 < 2e-16 ***
## k 37.03352 1.05271 35.179 < 2e-16 ***
## p 0.20173 0.01784 11.310 < 2e-16 ***
## s 2.64792 0.19044 13.904 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.3943 on 5177 degrees of freedom
##
## Number of iterations to convergence: 6
## Achieved convergence tolerance: 1.479e-06
## (3 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 5178 862.65
## 2 5177 803.67 1 58.977 379.91 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2c 56899.12
## 2 4 57255.49
## 3 5 56890.45
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.87572 0.14097 6.212 5.64e-10 ***
## alpha 0.82430 0.03969 20.766 < 2e-16 ***
## a 31.76130 2.62008 12.122 < 2e-16 ***
## b 122.35089 4.61204 26.529 < 2e-16 ***
## c 103.12833 4.07801 25.289 < 2e-16 ***
## d 1.27367 0.05860 21.735 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.394 on 5177 degrees of freedom
##
## Algorithm "port", convergence message: relative convergence (4)
## (3 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 596 80.866
## 2 595 71.900 1 8.9663 74.2 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 6115.535
## 2 2 6047.140
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau -0.07855 0.26298 -0.299 0.765
## alpha 0.89957 0.09748 9.228 < 2e-16 ***
## A 298.68647 41.73938 7.156 2.45e-12 ***
## k 95.82372 18.87728 5.076 5.16e-07 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.3476 on 595 degrees of freedom
##
## Number of iterations to convergence: 5
## Achieved convergence tolerance: 5.77e-06
## (3 observations deleted due to missingness)
## Error in nls(get(paste("f_", Mod.Sel1, "c", sep = "")), data = G_M223, :
## number of iterations exceeded maximum of 50
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 595 71.900
## 2 594 71.708 1 0.19119 1.5837 0.2087
## model AIC
## 1 2 6047.140
## 2 2a 6047.545
## 3 2b 6048.497
## 4 2c NA
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau -0.07855 0.26298 -0.299 0.765
## alpha 0.89957 0.09748 9.228 < 2e-16 ***
## A 298.68647 41.73938 7.156 2.45e-12 ***
## k 95.82372 18.87728 5.076 5.16e-07 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.3476 on 595 degrees of freedom
##
## Number of iterations to convergence: 5
## Achieved convergence tolerance: 5.77e-06
## (3 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 594 80.723
## 2 593 71.768 1 8.9544 73.988 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 6047.140
## 2 4 6118.473
## 3 5 6050.044
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau -0.07855 0.26298 -0.299 0.765
## alpha 0.89957 0.09748 9.228 < 2e-16 ***
## A 298.68647 41.73938 7.156 2.45e-12 ***
## k 95.82372 18.87728 5.076 5.16e-07 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.3476 on 595 degrees of freedom
##
## Number of iterations to convergence: 5
## Achieved convergence tolerance: 5.77e-06
## (3 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 676 152.77
## 2 675 142.27 1 10.498 49.809 4.216e-12 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 7025.553
## 2 2 6979.211
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.5697 0.5042 1.130 0.259
## alpha 0.8583 0.1139 7.536 1.57e-13 ***
## A 315.8268 62.8438 5.026 6.44e-07 ***
## k 147.8427 33.1149 4.465 9.40e-06 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4591 on 675 degrees of freedom
##
## Number of iterations to convergence: 5
## Achieved convergence tolerance: 4.062e-06
## (1 observation deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 675 142.27
## 2 674 142.26 1 0.0098593 0.0467 0.8290
## 3 674 142.26 0 0.0000000
## 4 673 142.26 1 0.0069224 0.0327 0.8564
## model AIC
## 1 2 6979.211
## 2 2a 6981.164
## 3 2b 6981.181
## 4 2c 6983.148
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.5697 0.5042 1.130 0.259
## alpha 0.8583 0.1139 7.536 1.57e-13 ***
## A 315.8268 62.8438 5.026 6.44e-07 ***
## k 147.8427 33.1149 4.465 9.40e-06 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4591 on 675 degrees of freedom
##
## Number of iterations to convergence: 5
## Achieved convergence tolerance: 4.062e-06
## (1 observation deleted due to missingness)
## Error in nls(f_5, data = G_M231, start = c(tau = tau.start, alpha = alpha.start, :
## Convergence failure: iteration limit reached without convergence (10)
## model AIC
## 1 2 6979.211
## 2 4 7029.273
## 3 5 NA
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.5697 0.5042 1.130 0.259
## alpha 0.8583 0.1139 7.536 1.57e-13 ***
## A 315.8268 62.8438 5.026 6.44e-07 ***
## k 147.8427 33.1149 4.465 9.40e-06 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4591 on 675 degrees of freedom
##
## Number of iterations to convergence: 5
## Achieved convergence tolerance: 4.062e-06
## (1 observation deleted due to missingness)
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 326 100.337
## 2 325 98.647 1 1.6905 5.5696 0.01887 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 3504.259
## 2 2 3500.669
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.5690 1.1894 1.319 0.18805
## alpha 0.6751 0.2661 2.537 0.01165 *
## A 85.0151 18.2427 4.660 4.61e-06 ***
## k 22.9171 7.8626 2.915 0.00381 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5509 on 325 degrees of freedom
##
## Number of iterations to convergence: 12
## Achieved convergence tolerance: 7.882e-06
## (1 observation deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 325 98.647
## 2 324 97.539 1 1.10741 3.6785 0.0560 .
## 3 324 93.712 0 0.00000
## 4 323 93.308 1 0.40411 1.3989 0.2378
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 3500.669
## 2 2a 3498.954
## 3 2b 3485.785
## 4 2c 3486.363
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.8945 1.2923 1.466 0.1436
## alpha 0.6580 0.2763 2.381 0.0178 *
## A 65.9248 13.3842 4.926 1.34e-06 ***
## k 30.8188 3.7464 8.226 4.75e-15 ***
## s 42.5340 187.7406 0.227 0.8209
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5378 on 324 degrees of freedom
##
## Number of iterations to convergence: 21
## Achieved convergence tolerance: 8.771e-06
## (1 observation deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 324 96.856
## 2 323 95.476 1 1.38 4.6688 0.03145 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2b 3485.785
## 2 4 3496.641
## 3 5 3493.920
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.8945 1.2923 1.466 0.1436
## alpha 0.6580 0.2763 2.381 0.0178 *
## A 65.9248 13.3842 4.926 1.34e-06 ***
## k 30.8188 3.7464 8.226 4.75e-15 ***
## s 42.5340 187.7406 0.227 0.8209
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5378 on 324 degrees of freedom
##
## Number of iterations to convergence: 21
## Achieved convergence tolerance: 8.771e-06
## (1 observation deleted due to missingness)
simple log-normal model: does not fit
log-normal phi model: does not fit
model can fit - but K is negative (only 19 observations) - model excluded
## [1] "cannot plot data with prediction"
## [1] "cannot plot data with prediction"
## [1] "cannot plot data with prediction"
## [1] "cannot plot data with prediction"
## [1] "cannot plot data with prediction"
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 299 124.35
## 2 298 113.37 1 10.974 28.846 1.581e-07 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 3013.807
## 2 2 2987.903
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.07662 0.66159 0.116 0.90788
## alpha 0.77829 0.12977 5.998 5.79e-09 ***
## A 99.54181 19.23214 5.176 4.18e-07 ***
## k 53.05797 17.65911 3.005 0.00289 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.6168 on 298 degrees of freedom
##
## Number of iterations to convergence: 5
## Achieved convergence tolerance: 6.942e-06
## (4 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 298 113.37
## 2 297 113.26 1 0.115630 0.3032 0.5823
## 3 297 113.00 0 0.000000
## 4 296 112.91 1 0.094618 0.2480 0.6188
## model AIC
## 1 2 2987.903
## 2 2a 2989.595
## 3 2b 2988.916
## 4 2c 2990.663
## Warning in `[<-.data.frame`(`*tmp*`, nls.param.df_bal$Code == "M334", , :
## provided 33 variables to replace 32 variables
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.07662 0.66159 0.116 0.90788
## alpha 0.77829 0.12977 5.998 5.79e-09 ***
## A 99.54181 19.23214 5.176 4.18e-07 ***
## k 53.05797 17.65911 3.005 0.00289 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.6168 on 298 degrees of freedom
##
## Number of iterations to convergence: 5
## Achieved convergence tolerance: 6.942e-06
## (4 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 297 124.07
## 2 296 113.29 1 10.775 28.151 2.203e-07 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 2987.903
## 2 4 3017.126
## 3 5 2991.690
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.07662 0.66159 0.116 0.90788
## alpha 0.77829 0.12977 5.998 5.79e-09 ***
## A 99.54181 19.23214 5.176 4.18e-07 ***
## k 53.05797 17.65911 3.005 0.00289 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.6168 on 298 degrees of freedom
##
## Number of iterations to convergence: 5
## Achieved convergence tolerance: 6.942e-06
## (4 observations deleted due to missingness)
## [1] "cannot plot data with prediction"
| Code | Ecoregion | Sel.Mod.2 | Sel.Mod.3 | Best.Mod |
|---|---|---|---|---|
| 211 | Northeastern Mixed Forest | 2 | 5 | 5 |
| 212 | Laurentian Mixed Forest | 2a | 5 | 5 |
| 221 | Eastern Broadleaf Forest | 2b | 5 | 5 |
| 222 | Midwest Broadleaf Forest | 2a | 5 | 5 |
| 223 | Central Interior Broadleaf Forest | 2c | 5 | 5 |
| 231 | Southeastern Mixed Forest | 2c | 5 | 5 |
| 232 | Outer Coastal Plain Mixed Forest | 2a | 5 | 5 |
| 234 | Lower Mississippi Riverine Forest | 2c | 5 | 5 |
| 242 | Pacific Lowland Mixed Forest | NA | NA | NA |
| 251 | Prairie Parkland (Temperate) | 2b | 5 | 5 |
| 255 | Prairie Parkland (Subtropical) | 2 | 5 | 5 |
| 261 | California Coastal Chaparral Forest and Shrub | NA | NA | NA |
| 262 | California Dry Steppe | NA | NA | NA |
| 263 | California Coastal Steppe - Mixed Forest and Redwood Forest | NA | NA | NA |
| 313 | Colorado Plateau Semi-Desert | NA | NA | NA |
| 315 | Southwest Plateau and Plains Dry Steppe and Shrub | NA | NA | NA |
| 321 | Chihuahuan Semi-Desert | NA | NA | NA |
| 322 | American Semidesert and Desert | NA | NA | NA |
| 331 | Great Plains/Palouse Dry Steppe | NA | NA | NA |
| 332 | Great Plains Steppe | 2 | 2 | 2 |
| 341 | Intermountain Semi-Desert and Desert | NA | NA | NA |
| 342 | Intermountain Semi-Desert | NA | NA | NA |
| 411 | Everglades | NA | NA | NA |
| M211 | Adirondack-New England Mixed forest - Coniferous Forest - Alpine Meadow | 2a | 5 | 5 |
| M221 | Central Appalachian Broadleaf Forest - Coniferous Forest - Meadow | 2c | 5 | 5 |
| M223 | Ozark Broadleaf Forest Meadow | 2 | 2 | 2 |
| M231 | Ouachita Mixed Forest | 2 | 2 | 2 |
| M242 | Cascade Mixed Forest | NA | NA | NA |
| M261 | Sierran Steppe - Mixed Forest - Coniferous Forest - Alpine Meadow | 2b | 2b | 2b |
| M262 | California Coastal Range Coniferous Forest - Open Woodland - Shrub - Meadow | NA | NA | NA |
| M313 | Arizona-New Mexico Mountains Semi-Desert - Open Woodland - Coniferous Forest - Alpine Meadow | NA | NA | NA |
| M331 | Southern Rocky Mountain Steppe - Open Woodland - Coniferous Forest - Alpine Meadow | NA | NA | NA |
| M332 | Middle Rocky Mountain Steppe - Coniferous Forest - Alpine Meadow | NA | NA | NA |
| M333 | Northern Rocky Mountain Steppe - Coniferous Forest - Alpine Meadow | NA | NA | NA |
| M334 | Black Hills Coniferous Forest | 2 | 2 | 2 |
| M341 | Nevada-Utah Mountains Semi-Desert - Coniferous Forest - Alpine Meadow | NA | NA | NA |
| Code | Ecoregion | region | n.obs | n.plots | tau | tau.variance | tau.2.5 | tau.97.5 | alpha | alpha.variance | alpha.2.5 | alpha.97.5 | A | A.2.5 | A.97.5 | k | k.2.5 | k.97.5 | a | a.2.5 | a.97.5 | b | b.se | b.2.5 | b.97.5 | c | c.2.5 | c.97.5 | d | d.2.5 | d.97.5 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 211 | Northeastern Mixed Forest | east | 4838 | 2419 | 0.5811807 | 0.0315130 | 0.2331627 | 0.9291986 | 0.8407842 | 0.0011638 | 0.7739056 | 0.9076629 | 473.76484 | 401.22103 | 546.30865 | 205.69316 | 169.97959 | 241.40674 | 33.400705 | 29.493797 | 37.30761 | 118.88478 | NA | 107.26593 | 130.50363 | 118.37512 | 106.88772 | 129.86252 | 1.016528 | 0.9139327 | 1.119122 |
| 212 | Laurentian Mixed Forest | east | 12962 | 6481 | 1.1905141 | 0.0185145 | 0.9238008 | 1.4572274 | 0.6733287 | 0.0007299 | 0.6203705 | 0.7262870 | 192.14309 | 173.82377 | 210.46242 | 108.46141 | 92.82129 | 124.10153 | 13.385905 | 12.379552 | 14.39226 | 75.80154 | NA | 71.16321 | 80.43987 | 121.05709 | 110.36363 | 131.75055 | 1.415947 | 1.3319180 | 1.499977 |
| 221 | Eastern Broadleaf Forest | east | 5446 | 2723 | 0.2422230 | 0.0169479 | -0.0129899 | 0.4974359 | 0.8082238 | 0.0009532 | 0.7476984 | 0.8687492 | 293.53650 | 251.49470 | 335.57830 | 64.20469 | 51.62403 | 76.78536 | 21.082556 | 15.831554 | 26.33356 | 179.51639 | NA | 159.02822 | 200.00457 | 156.83930 | 126.03502 | 187.64358 | 1.518786 | 1.3427586 | 1.694814 |
| 222 | Midwest Broadleaf Forest | east | 3552 | 1776 | 1.2114537 | 0.0651252 | 0.7111070 | 1.7118003 | 0.7688102 | 0.0029588 | 0.6621620 | 0.8754583 | 342.27169 | 256.80774 | 427.73563 | 170.04542 | 113.41537 | 226.67548 | 14.018976 | 11.290384 | 16.74757 | 101.92211 | NA | 90.10999 | 113.73424 | 115.56990 | 99.46147 | 131.67834 | 1.197007 | 1.0662985 | 1.327715 |
| 223 | Central Interior Broadleaf Forest | east | 6388 | 3194 | 0.8461273 | 0.0146625 | 0.6087525 | 1.0835021 | 0.7112156 | 0.0011445 | 0.6448959 | 0.7775352 | 127.11413 | 118.06312 | 136.16514 | 37.02397 | 34.56160 | 39.48633 | 17.800944 | 13.742432 | 21.85946 | 98.51432 | NA | 90.85761 | 106.17104 | 113.18502 | 99.49412 | 126.87591 | 1.445696 | 1.2995960 | 1.591796 |
| 231 | Southeastern Mixed Forest | east | 7790 | 3895 | 1.4759858 | 0.0331065 | 1.1193112 | 1.8326605 | 0.6508476 | 0.0005086 | 0.6066399 | 0.6950552 | 144.97930 | 132.73182 | 157.22678 | 30.16077 | 27.76581 | 32.55574 | 16.127853 | 14.505323 | 17.75038 | 120.36292 | NA | 109.42373 | 131.30210 | 112.83900 | 94.67141 | 131.00659 | 1.633502 | 1.5036347 | 1.763370 |
| 232 | Outer Coastal Plain Mixed Forest | east | 7940 | 3970 | 1.0085159 | 0.0364082 | 0.6344793 | 1.3825526 | 0.7104306 | 0.0004570 | 0.6685254 | 0.7523358 | 408.34768 | 328.15225 | 488.54310 | 159.56714 | 117.80214 | 201.33213 | 20.518658 | 18.511922 | 22.52539 | 121.27486 | NA | 108.92335 | 133.62637 | 111.11452 | 92.96303 | 129.26601 | 1.502835 | 1.3689187 | 1.636751 |
| 234 | Lower Mississippi Riverine Forest | east | 830 | 415 | 0.4023597 | 0.1697654 | -0.4063858 | 1.2111052 | 0.7883253 | 0.0049140 | 0.6507290 | 0.9259216 | 395.47217 | 23.90610 | 767.03825 | 99.96670 | -38.80593 | 238.73932 | 0.000000 | -11.644423 | 11.64442 | 340.82383 | NA | -60.32947 | 741.97714 | 675.40707 | -1170.66705 | 2521.48119 | 2.554711 | 1.2035154 | 3.905907 |
| 242 | Pacific Lowland Mixed Forest | west | 0 | 0 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 251 | Prairie Parkland (Temperate) | east | 1392 | 696 | 0.3416233 | 0.0635084 | -0.1527370 | 0.8359836 | 0.7012428 | 0.0088949 | 0.5162318 | 0.8862538 | 142.90894 | 116.45330 | 169.36458 | 39.02880 | 31.20558 | 46.85201 | 24.394348 | 16.444251 | 32.34445 | 93.28179 | NA | 78.72633 | 107.83724 | 102.73710 | 85.23504 | 120.23915 | 1.135955 | 0.9228090 | 1.349101 |
| 255 | Prairie Parkland (Subtropical) | east | 444 | 222 | -0.2439613 | 0.3026388 | -1.3251902 | 0.8372676 | 0.5283394 | 0.0248661 | 0.2184128 | 0.8382661 | 320.37518 | 110.66455 | 530.08580 | 143.22594 | 39.60855 | 246.84333 | 9.564518 | 5.735244 | 13.39379 | 88.51384 | NA | 63.99654 | 113.03113 | 55.99135 | 44.35208 | 67.63061 | 1.039193 | 0.8154724 | 1.262914 |
| 261 | California Coastal Chaparral Forest and Shrub | west | 0 | 0 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 262 | California Dry Steppe | west | 0 | 0 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 263 | California Coastal Steppe - Mixed Forest and Redwood Forest | west | 4 | 2 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 313 | Colorado Plateau Semi-Desert | west | 0 | 0 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 315 | Southwest Plateau and Plains Dry Steppe and Shrub | west | 0 | 0 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 321 | Chihuahuan Semi-Desert | west | 0 | 0 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 322 | American Semidesert and Desert | west | 0 | 0 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 331 | Great Plains/Palouse Dry Steppe | west | 118 | 59 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 332 | Great Plains Steppe | west | 154 | 77 | 3.0054586 | 17.0455900 | -5.1532211 | 11.1641382 | 0.9713883 | 0.0886175 | 0.3831227 | 1.5596540 | 231.92537 | -216.63238 | 680.48311 | 219.82421 | -175.36661 | 615.01504 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 341 | Intermountain Semi-Desert and Desert | west | 4 | 2 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 342 | Intermountain Semi-Desert | west | 2 | 1 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 411 | Everglades | east | 66 | 33 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| M211 | Adirondack-New England Mixed forest - Coniferous Forest - Alpine Meadow | east | 5108 | 2554 | 0.5302798 | 0.0245169 | 0.2233184 | 0.8372412 | 0.8138994 | 0.0007135 | 0.7615330 | 0.8662658 | 324.77328 | 285.05180 | 364.49475 | 102.66625 | 81.59293 | 123.73957 | 14.565033 | 8.704894 | 20.42517 | 155.65092 | NA | 135.84666 | 175.45518 | 179.79825 | 140.99259 | 218.60392 | 1.586471 | 1.3826673 | 1.790274 |
| M221 | Central Appalachian Broadleaf Forest - Coniferous Forest - Meadow | east | 5186 | 2593 | 0.8757158 | 0.0198734 | 0.5993491 | 1.1520824 | 0.8243016 | 0.0015756 | 0.7464844 | 0.9021187 | 162.45200 | 151.44589 | 173.45812 | 37.03352 | 34.96976 | 39.09727 | 31.761299 | 26.624837 | 36.89776 | 122.35089 | NA | 113.30935 | 131.39243 | 103.12833 | 95.13370 | 111.12296 | 1.273667 | 1.1587869 | 1.388547 |
| M223 | Ozark Broadleaf Forest Meadow | east | 602 | 301 | -0.0785477 | 0.0691598 | -0.5950347 | 0.4379392 | 0.8995732 | 0.0095032 | 0.7081176 | 1.0910287 | 298.68647 | 216.71203 | 380.66091 | 95.82372 | 58.74951 | 132.89792 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| M231 | Ouachita Mixed Forest | east | 680 | 340 | 0.5697483 | 0.2542625 | -0.4203280 | 1.5598245 | 0.8582543 | 0.0129720 | 0.6346241 | 1.0818845 | 315.82677 | 192.43391 | 439.21963 | 147.84274 | 82.82219 | 212.86329 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| M242 | Cascade Mixed Forest | west | 34 | 17 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| M261 | Sierran Steppe - Mixed Forest - Coniferous Forest - Alpine Meadow | west | 330 | 165 | 1.8945480 | 1.6701111 | -0.6478649 | 4.4369609 | 0.6580466 | 0.0763610 | 0.1144095 | 1.2016837 | 65.92482 | 39.59387 | 92.25577 | 30.81882 | 23.44844 | 38.18920 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| M262 | California Coastal Range Coniferous Forest - Open Woodland - Shrub - Meadow | west | 8 | 4 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| M313 | Arizona-New Mexico Mountains Semi-Desert - Open Woodland - Coniferous Forest - Alpine Meadow | west | 0 | 0 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| M331 | Southern Rocky Mountain Steppe - Open Woodland - Coniferous Forest - Alpine Meadow | west | 0 | 0 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| M332 | Middle Rocky Mountain Steppe - Coniferous Forest - Alpine Meadow | west | 20 | 10 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| M333 | Northern Rocky Mountain Steppe - Coniferous Forest - Alpine Meadow | west | 22 | 11 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| M334 | Black Hills Coniferous Forest | west | 306 | 153 | 0.0766154 | 0.4376976 | -1.2253593 | 1.3785901 | 0.7782879 | 0.0168399 | 0.5229091 | 1.0336667 | 99.54181 | 61.69380 | 137.38982 | 53.05797 | 18.30561 | 87.81033 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| M341 | Nevada-Utah Mountains Semi-Desert - Coniferous Forest - Alpine Meadow | west | 0 | 0 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
## OGR data source with driver: ESRI Shapefile
## Source: "C:\Users\hogan.jaaron\Dropbox\FIA_R\Mapping\S_USA.EcoMapProvinces\S_USA.EcoMapProvinces.shp", layer: "S_USA.EcoMapProvinces"
## with 37 features
## It has 17 fields
## Integer64 fields read as strings: PROVINCE_ PROVINCE_I
## Warning: Removed 19 rows containing missing values (`geom_point()`).
## Warning: Removed 19 rows containing missing values (`geom_point()`).
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 4767 1100.1
## 2 4766 1043.0 1 57.113 260.99 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 51935.28
## 2 2 51682.97
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.55012 0.29122 5.323 1.07e-07 ***
## alpha 1.10900 0.06402 17.322 < 2e-16 ***
## A 358.60119 28.52883 12.570 < 2e-16 ***
## k 179.16216 14.20596 12.612 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4678 on 4766 degrees of freedom
##
## Number of iterations to convergence: 7
## Achieved convergence tolerance: 2.058e-06
## (2 observations deleted due to missingness)
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 4766 1043
## 2 4765 1041 1 1.9622 8.982 0.002741 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 51682.97
## 2 2a 51675.99
## 3 2b 51658.42
## 4 2c NA
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.45638 0.28154 5.173 2.40e-07 ***
## alpha 1.12704 0.06463 17.440 < 2e-16 ***
## A 214.23259 19.48572 10.994 < 2e-16 ***
## k 72.38593 8.79876 8.227 2.46e-16 ***
## s 1.30993 0.06714 19.511 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4665 on 4765 degrees of freedom
##
## Number of iterations to convergence: 8
## Achieved convergence tolerance: 8.183e-06
## (2 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 4765 1083.1
## 2 4764 1027.3 1 55.751 258.53 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2b 51658.42
## 2 4 51865.15
## 3 5 51615.07
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.58266 0.29313 5.399 7.02e-08 ***
## alpha 1.10888 0.06432 17.240 < 2e-16 ***
## a 23.88706 2.03910 11.714 < 2e-16 ***
## b 108.49559 7.02793 15.438 < 2e-16 ***
## c 132.60613 8.82813 15.021 < 2e-16 ***
## d 1.19405 0.06821 17.506 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4644 on 4764 degrees of freedom
##
## Algorithm "port", convergence message: relative convergence (4)
## (2 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 17134 6198.8
## 2 17133 6050.3 1 148.5 420.51 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 178725.1
## 2 2 178311.6
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.10028 0.13197 8.337 <2e-16 ***
## alpha 0.81927 0.03764 21.768 <2e-16 ***
## A 152.07486 4.72755 32.168 <2e-16 ***
## k 73.19820 2.56503 28.537 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5943 on 17133 degrees of freedom
##
## Number of iterations to convergence: 5
## Achieved convergence tolerance: 8.52e-06
## (13 observations deleted due to missingness)
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 17133 6050.3
## 2 17132 6047.1 1 3.2808 9.2947 0.002302 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 178311.6
## 2 2a 178304.3
## 3 2b 178234.1
## 4 2c NA
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.16176 0.13429 8.651 <2e-16 ***
## alpha 0.83619 0.03798 22.015 <2e-16 ***
## A 112.92651 4.28367 26.362 <2e-16 ***
## k 41.69107 2.16303 19.274 <2e-16 ***
## s 1.30455 0.03898 33.467 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5929 on 17132 degrees of freedom
##
## Number of iterations to convergence: 7
## Achieved convergence tolerance: 9.046e-06
## (13 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 17132 6118.8
## 2 17131 5964.1 1 154.7 444.35 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2b 178234.1
## 2 4 178506.3
## 3 5 178069.5
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.21466 0.13596 8.934 <2e-16 ***
## alpha 0.82392 0.03678 22.402 <2e-16 ***
## a 13.98018 0.80587 17.348 <2e-16 ***
## b 71.24841 2.09216 34.055 <2e-16 ***
## c 116.92499 4.25601 27.473 <2e-16 ***
## d 1.41175 0.04295 32.869 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.59 on 17131 degrees of freedom
##
## Algorithm "port", convergence message: relative convergence (4)
## (13 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 5798 981.08
## 2 5797 926.82 1 54.253 339.34 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 63460.21
## 2 2 63132.20
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.51229 0.15326 3.343 0.000835 ***
## alpha 0.88535 0.04543 19.490 < 2e-16 ***
## A 508.99853 32.76656 15.534 < 2e-16 ***
## k 174.84347 12.40713 14.092 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.3998 on 5797 degrees of freedom
##
## Number of iterations to convergence: 7
## Achieved convergence tolerance: 5.849e-06
## (1 observation deleted due to missingness)
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 5797 926.82
## 2 5796 919.97 1 6.8551 43.189 5.405e-11 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 63132.20
## 2 2a 63091.14
## 3 2b 58973.93
## 4 2c NA
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.23383 0.13009 1.797 0.0723 .
## alpha 0.81624 0.03142 25.975 <2e-16 ***
## A 293.53650 21.44552 13.688 <2e-16 ***
## k 64.20469 6.41740 10.005 <2e-16 ***
## s 1.31952 0.05962 22.132 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.3959 on 5440 degrees of freedom
##
## Number of iterations to convergence: 7
## Achieved convergence tolerance: 2.275e-06
## (1 observation deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 5796 965.51
## 2 5795 910.77 1 54.739 348.29 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2b 58973.93
## 2 4 63371.45
## 3 5 63034.87
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.23383 0.13009 1.797 0.0723 .
## alpha 0.81624 0.03142 25.975 <2e-16 ***
## A 293.53650 21.44552 13.688 <2e-16 ***
## k 64.20469 6.41740 10.005 <2e-16 ***
## s 1.31952 0.05962 22.132 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.3959 on 5440 degrees of freedom
##
## Number of iterations to convergence: 7
## Achieved convergence tolerance: 2.275e-06
## (1 observation deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 4359 1525.4
## 2 4358 1473.0 1 52.402 155.04 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 46875.15
## 2 2 46724.67
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.79200 0.23894 3.315 0.000925 ***
## alpha 0.84551 0.06349 13.318 < 2e-16 ***
## A 630.56526 94.45976 6.675 2.77e-11 ***
## k 370.30471 59.25802 6.249 4.53e-10 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5814 on 4358 degrees of freedom
##
## Number of iterations to convergence: 8
## Achieved convergence tolerance: 7.781e-06
## (2 observations deleted due to missingness)
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 4358 1473.0
## 2 4357 1460.9 1 12.03 35.879 2.269e-09 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 46724.67
## 2 2a 46690.90
## 3 2b NA
## 4 2c NA
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.865508 0.243990 3.547 0.000393 ***
## alpha 0.877804 0.063130 13.905 < 2e-16 ***
## A 382.711524 46.064067 8.308 < 2e-16 ***
## k 181.086971 27.994519 6.469 1.1e-10 ***
## p -0.021832 0.006205 -3.518 0.000439 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5791 on 4357 degrees of freedom
##
## Number of iterations to convergence: 7
## Achieved convergence tolerance: 8.562e-06
## (2 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 4357 1488.5
## 2 4356 1424.2 1 64.342 196.8 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2a 46690.90
## 2 4 46772.49
## 3 5 46581.74
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.00263 0.25225 3.975 7.16e-05 ***
## alpha 0.92051 0.06062 15.184 < 2e-16 ***
## a 12.37804 1.27551 9.704 < 2e-16 ***
## b 109.05738 6.07768 17.944 < 2e-16 ***
## c 113.94899 6.97245 16.343 < 2e-16 ***
## d 1.18472 0.05679 20.861 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5718 on 4356 degrees of freedom
##
## Algorithm "port", convergence message: relative convergence (4)
## (2 observations deleted due to missingness)
##
## ------
## Shapiro-Wilk normality test
##
## data: stdres
## W = 0.96113, p-value < 2.2e-16
##
##
## ------
##
## Runs Test
##
## data: as.factor(run)
## Standard Normal = -31.622, p-value < 2.2e-16
## alternative hypothesis: two.sided
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 7425 1362.9
## 2 7424 1316.5 1 46.482 262.13 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 77578.06
## 2 2 77322.31
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.04382 0.14030 7.44 1.12e-13 ***
## alpha 0.79468 0.04635 17.15 < 2e-16 ***
## A 229.08813 9.44538 24.25 < 2e-16 ***
## k 88.54009 4.45866 19.86 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4211 on 7424 degrees of freedom
##
## Number of iterations to convergence: 7
## Achieved convergence tolerance: 1.902e-06
## (3 observations deleted due to missingness)
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 7424 1316.5
## 2 7423 1312.4 1 4.078 23.0681 1.594e-06 ***
## 3 7423 1296.2 0 0.000
## 4 6380 1128.1 1043 168.074 0.9114 0.9733
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 77322.31
## 2 2a 77301.26
## 3 2b 77208.80
## 4 2c 66305.16
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.85759 0.12173 7.045 2.05e-12 ***
## alpha 0.71209 0.03381 21.062 < 2e-16 ***
## A 127.11413 4.61707 27.531 < 2e-16 ***
## k 37.02397 1.25609 29.475 < 2e-16 ***
## p 0.14094 0.01883 7.486 8.09e-14 ***
## s 2.17271 0.15279 14.220 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4205 on 6380 degrees of freedom
##
## Number of iterations to convergence: 5
## Achieved convergence tolerance: 4.063e-06
## (2 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 7423 1331.0
## 2 7422 1284.9 1 46.068 266.1 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2c 66305.16
## 2 4 77405.77
## 3 5 77146.11
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.85759 0.12173 7.045 2.05e-12 ***
## alpha 0.71209 0.03381 21.062 < 2e-16 ***
## A 127.11413 4.61707 27.531 < 2e-16 ***
## k 37.02397 1.25609 29.475 < 2e-16 ***
## p 0.14094 0.01883 7.486 8.09e-14 ***
## s 2.17271 0.15279 14.220 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4205 on 6380 degrees of freedom
##
## Number of iterations to convergence: 5
## Achieved convergence tolerance: 4.063e-06
## (2 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 9262 2333.1
## 2 9261 2269.7 1 63.327 258.39 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 101389.5
## 2 2 101136.6
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.58614 0.16735 9.478 <2e-16 ***
## alpha 0.71207 0.04222 16.865 <2e-16 ***
## A 286.74176 11.66804 24.575 <2e-16 ***
## k 85.62977 3.37714 25.356 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4951 on 9261 degrees of freedom
##
## Number of iterations to convergence: 6
## Achieved convergence tolerance: 1.425e-06
## (2 observations deleted due to missingness)
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 9261 2269.7
## 2 9260 2263.6 1 6.179 25.2788 5.054e-07 ***
## 3 7771 2566.4 1489 -302.856 -0.6159 1
## 4 7770 2496.6 1 69.841 217.3639 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 101136.58
## 2 2a 101113.33
## 3 2b 85459.40
## 4 2c 85246.85
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.506e+00 1.836e-01 8.201 2.76e-16 ***
## alpha 6.512e-01 2.258e-02 28.840 < 2e-16 ***
## A 1.450e+02 6.248e+00 23.205 < 2e-16 ***
## k 3.016e+01 1.222e+00 24.686 < 2e-16 ***
## p 1.076e-01 7.674e-03 14.024 < 2e-16 ***
## s 1.963e+00 9.744e-02 20.143 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5668 on 7770 degrees of freedom
##
## Number of iterations to convergence: 12
## Achieved convergence tolerance: 5.033e-06
## (14 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 9260 2303.2
## 2 9259 2236.8 1 66.418 274.94 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2c 85246.85
## 2 4 101274.10
## 3 5 101004.99
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.506e+00 1.836e-01 8.201 2.76e-16 ***
## alpha 6.512e-01 2.258e-02 28.840 < 2e-16 ***
## A 1.450e+02 6.248e+00 23.205 < 2e-16 ***
## k 3.016e+01 1.222e+00 24.686 < 2e-16 ***
## p 1.076e-01 7.674e-03 14.024 < 2e-16 ***
## s 1.963e+00 9.744e-02 20.143 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5668 on 7770 degrees of freedom
##
## Number of iterations to convergence: 12
## Achieved convergence tolerance: 5.033e-06
## (14 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 9495 3706.9
## 2 9494 3585.7 1 121.15 320.78 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 106844.7
## 2 2 106531.1
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.12108 0.18182 6.166 7.29e-10 ***
## alpha 0.74475 0.03909 19.052 < 2e-16 ***
## A 303.18915 15.39059 19.700 < 2e-16 ***
## k 91.43033 4.53472 20.162 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.6146 on 9494 degrees of freedom
##
## Number of iterations to convergence: 5
## Achieved convergence tolerance: 6.341e-06
## (8 observations deleted due to missingness)
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 9494 3585.7
## 2 9493 3580.8 1 4.8991 12.988 0.0003151 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 106531.09
## 2 2a 106520.11
## 3 2b 88567.74
## 4 2c NA
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.83554 0.18434 4.533 5.91e-06 ***
## alpha 0.58429 0.02782 21.001 < 2e-16 ***
## A 700.20538 245.84610 2.848 0.00441 **
## k 453.42710 284.16976 1.596 0.11061
## s 0.77872 0.03365 23.145 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.7037 on 7903 degrees of freedom
##
## Number of iterations to convergence: 12
## Achieved convergence tolerance: 2.783e-06
## (32 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 9493 3686.3
## 2 9492 3559.9 1 126.38 336.99 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2b 88567.74
## 2 4 106795.73
## 3 5 106466.38
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.83554 0.18434 4.533 5.91e-06 ***
## alpha 0.58429 0.02782 21.001 < 2e-16 ***
## A 700.20538 245.84610 2.848 0.00441 **
## k 453.42710 284.16976 1.596 0.11061
## s 0.77872 0.03365 23.145 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.7037 on 7903 degrees of freedom
##
## Number of iterations to convergence: 12
## Achieved convergence tolerance: 2.783e-06
## (32 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 1083 356.24
## 2 1082 339.64 1 16.601 52.886 6.77e-13 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 12290.56
## 2 2 12240.74
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.9025 0.6136 1.471 0.1417
## alpha 0.8625 0.1098 7.854 9.68e-15 ***
## A 1162.2886 522.4216 2.225 0.0263 *
## k 500.8135 237.4249 2.109 0.0351 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5603 on 1082 degrees of freedom
##
## Number of iterations to convergence: 7
## Achieved convergence tolerance: 2.518e-06
## (9 observations deleted due to missingness)
## Error in nls(get(paste("f_", Mod.Sel1, "a", sep = "")), data = G_234, :
## step factor 0.000488281 reduced below 'minFactor' of 0.000976562
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 1082 339.64
## 2 823 192.48 259 147.16 2.4293 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 12240.738
## 2 2a NA
## 3 2b NA
## 4 2c 8996.781
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.39536 0.41113 0.962 0.337
## alpha 0.78925 0.07021 11.242 < 2e-16 ***
## A 395.47217 189.29921 2.089 0.037 *
## k 99.96670 70.69953 1.414 0.158
## p -0.01220 0.01423 -0.857 0.392
## s 1.13109 0.24438 4.628 4.28e-06 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4836 on 823 degrees of freedom
##
## Number of iterations to convergence: 40
## Achieved convergence tolerance: 3.1e-06
## (1 observation deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 1081 355.31
## 2 1080 337.44 1 17.869 57.192 8.426e-14 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2c 8996.781
## 2 4 12291.729
## 3 5 12237.690
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.39536 0.41113 0.962 0.337
## alpha 0.78925 0.07021 11.242 < 2e-16 ***
## A 395.47217 189.29921 2.089 0.037 *
## k 99.96670 70.69953 1.414 0.158
## p -0.01220 0.01423 -0.857 0.392
## s 1.13109 0.24438 4.628 4.28e-06 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4836 on 823 degrees of freedom
##
## Number of iterations to convergence: 40
## Achieved convergence tolerance: 3.1e-06
## (1 observation deleted due to missingness)
##
## ------
## Shapiro-Wilk normality test
##
## data: stdres
## W = 0.95371, p-value = 1.756e-15
##
##
## ------
##
## Runs Test
##
## data: as.factor(run)
## Standard Normal = -10.597, p-value < 2.2e-16
## alternative hypothesis: two.sided
## [1] "cannot plot residuals"
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 1884 530.49
## 2 1883 520.40 1 10.091 36.511 1.825e-09 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 19835.95
## 2 2 19801.71
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.18712 0.25641 0.730 0.466
## alpha 0.63418 0.09964 6.364 2.45e-10 ***
## A 298.73479 32.98617 9.056 < 2e-16 ***
## k 126.75170 16.66267 7.607 4.41e-14 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5257 on 1883 degrees of freedom
##
## Number of iterations to convergence: 7
## Achieved convergence tolerance: 2.294e-06
## (3 observations deleted due to missingness)
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 1883 520.40
## 2 1882 520.18 1 0.21527 0.7788 0.3776
## model AIC
## 1 2 19801.71
## 2 2a 19802.93
## 3 2b 19791.60
## 4 2c NA
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.2497 0.2619 0.954 0.34
## alpha 0.6703 0.1001 6.699 2.76e-11 ***
## A 167.9953 18.6424 9.011 < 2e-16 ***
## k 47.7477 6.8077 7.014 3.22e-12 ***
## s 1.4586 0.1242 11.742 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5242 on 1882 degrees of freedom
##
## Number of iterations to convergence: 9
## Achieved convergence tolerance: 7.144e-06
## (3 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 1882 514.77
## 2 1881 503.84 1 10.929 40.8 2.119e-10 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2b 19791.60
## 2 4 19783.19
## 3 5 19744.70
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.20201 0.25282 0.799 0.424
## alpha 0.65588 0.09727 6.743 2.06e-11 ***
## a 21.61048 2.71796 7.951 3.16e-15 ***
## b 100.87424 7.44340 13.552 < 2e-16 ***
## c 107.78366 10.02239 10.754 < 2e-16 ***
## d 1.20475 0.09772 12.328 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5176 on 1881 degrees of freedom
##
## Algorithm "port", convergence message: relative convergence (4)
## (3 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 633 274.88
## 2 632 263.69 1 11.189 26.816 3.012e-07 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 6532.531
## 2 2 6508.102
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.1250 0.5533 0.226 0.821390
## alpha 0.7041 0.1238 5.687 1.98e-08 ***
## A 296.3420 70.9096 4.179 3.34e-05 ***
## k 128.7173 33.4109 3.853 0.000129 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.6459 on 632 degrees of freedom
##
## Number of iterations to convergence: 8
## Achieved convergence tolerance: 3.502e-06
## (2 observations deleted due to missingness)
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 632 263.69
## 2 631 262.60 1 1.0917 2.6232 0.1058
## model AIC
## 1 2 6508.102
## 2 2a 6507.463
## 3 2b NA
## 4 2c NA
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.08586 0.54208 0.158 0.87420
## alpha 0.73491 0.12324 5.963 4.11e-09 ***
## A 242.91208 59.23770 4.101 4.66e-05 ***
## k 88.31753 30.37767 2.907 0.00377 **
## p -0.01872 0.01659 -1.128 0.25970
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.6451 on 631 degrees of freedom
##
## Number of iterations to convergence: 13
## Achieved convergence tolerance: 6.185e-06
## (2 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 631 262.85
## 2 630 247.58 1 15.265 38.844 8.404e-10 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2a 6507.463
## 2 4 6508.057
## 3 5 6472.004
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau -0.08353 0.48167 -0.173 0.862
## alpha 0.79872 0.11384 7.016 5.91e-12 ***
## a 11.01478 2.30217 4.785 2.14e-06 ***
## b 88.30672 10.55299 8.368 3.79e-16 ***
## c 58.54755 5.05845 11.574 < 2e-16 ***
## d 1.06975 0.09500 11.261 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.6269 on 630 degrees of freedom
##
## Algorithm "port", convergence message: relative convergence (4)
## (2 observations deleted due to missingness)
add p model: does not fit
add s model: does not fit
add s+p model: does not fit
unable to fit model (only 64 observations)
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
add p model: does not fit
add s model: does not fit
add s+p model: does not fit
unable to fit model (0 observations)
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## Error in nls(f_1, data = G_322, start = c(tau = tau.start, A = A.start, :
## step factor 0.000488281 reduced below 'minFactor' of 0.000976562
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## model AIC
## 1 1 NA
## 2 2 NA
## Warning in min(AIC1_322$AIC, na.rm = T): no non-missing arguments to min;
## returning Inf
## Error in h(simpleError(msg, call)) :
## error in evaluating the argument 'object' in selecting a method for function 'summary': object 'nls_322.' not found
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 216 145.31
## 2 215 143.71 1 1.5947 2.3858 0.1239
## model AIC
## 1 1 2335.903
## 2 2 2335.487
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.5363 1.8691 0.287 0.774
## alpha 0.5955 0.3508 1.698 0.091 .
## A 1637.2126 6225.7800 0.263 0.793
## k 1408.2568 5450.4162 0.258 0.796
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.8176 on 215 degrees of freedom
##
## Number of iterations to convergence: 8
## Achieved convergence tolerance: 2.297e-06
## (2 observations deleted due to missingness)
## Error in nls(get(paste("f_", Mod.Sel1, "a", sep = "")), data = G_332, :
## step factor 0.000488281 reduced below 'minFactor' of 0.000976562
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 215 143.711
## 2 147 83.108 68 60.603 1.5764 0.01166 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 2335.487
## 2 2a NA
## 3 2b 1636.402
## 4 2c NA
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 2.8328 3.9193 0.723 0.47097
## alpha 0.9760 0.3039 3.211 0.00162 **
## A 85.6246 59.9298 1.429 0.15520
## k 44.9763 23.3407 1.927 0.05592 .
## s 1.7214 0.6592 2.611 0.00995 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.7519 on 147 degrees of freedom
##
## Number of iterations to convergence: 33
## Achieved convergence tolerance: 8.579e-06
## (2 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 214 137.79
## 2 213 135.66 1 2.1317 3.3469 0.06873 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2b 1636.402
## 2 4 2328.276
## 3 5 2326.862
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 2.8328 3.9193 0.723 0.47097
## alpha 0.9760 0.3039 3.211 0.00162 **
## A 85.6246 59.9298 1.429 0.15520
## k 44.9763 23.3407 1.927 0.05592 .
## s 1.7214 0.6592 2.611 0.00995 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.7519 on 147 degrees of freedom
##
## Number of iterations to convergence: 33
## Achieved convergence tolerance: 8.579e-06
## (2 observations deleted due to missingness)
simple log-normal model: does not fit
log-normal phi model: does not fit
model not fitted because only 62 observations
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 4834 775.91
## 2 4833 746.70 1 29.21 189.06 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 51220.60
## 2 2 51036.99
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.81073 0.28200 6.421 1.48e-10 ***
## alpha 0.87838 0.06026 14.577 < 2e-16 ***
## A 332.52443 23.77680 13.985 < 2e-16 ***
## k 173.44970 11.22937 15.446 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.3931 on 4833 degrees of freedom
##
## Number of iterations to convergence: 6
## Achieved convergence tolerance: 4.621e-06
## (1 observation deleted due to missingness)
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 4833 746.70
## 2 4832 738.81 1 7.8916 51.613 7.79e-13 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 51036.99
## 2 2a 50987.60
## 3 2b 50965.66
## 4 2c NA
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.52010 0.25515 5.958 2.74e-09 ***
## alpha 0.91818 0.06060 15.152 < 2e-16 ***
## A 181.37233 11.65322 15.564 < 2e-16 ***
## k 57.00052 4.15035 13.734 < 2e-16 ***
## s 1.46437 0.06059 24.167 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.3901 on 4832 degrees of freedom
##
## Number of iterations to convergence: 7
## Achieved convergence tolerance: 6.556e-06
## (1 observation deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 4832 765.27
## 2 4831 733.83 1 31.442 206.99 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2b 50965.66
## 2 4 51157.82
## 3 5 50956.90
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.56244 0.25921 6.028 1.79e-09 ***
## alpha 0.92196 0.06021 15.312 < 2e-16 ***
## a 12.85464 2.78039 4.623 3.87e-06 ***
## b 125.02105 8.62168 14.501 < 2e-16 ***
## c 166.06497 16.00981 10.373 < 2e-16 ***
## d 1.54602 0.09924 15.579 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.3897 on 4831 degrees of freedom
##
## Algorithm "port", convergence message: relative convergence (4)
## (1 observation deleted due to missingness)
##
## ------
## Shapiro-Wilk normality test
##
## data: stdres
## W = 0.98658, p-value < 2.2e-16
##
##
## ------
##
## Runs Test
##
## data: as.factor(run)
## Standard Normal = -33.203, p-value < 2.2e-16
## alternative hypothesis: two.sided
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 7196 1178.4
## 2 7195 1129.9 1 48.537 309.08 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 79716.90
## 2 2 79416.11
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.9250 0.1344 6.881 6.43e-12 ***
## alpha 0.9052 0.0490 18.474 < 2e-16 ***
## A 258.3782 8.9769 28.783 < 2e-16 ***
## k 63.0567 2.7132 23.241 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.3963 on 7195 degrees of freedom
##
## Number of iterations to convergence: 7
## Achieved convergence tolerance: 4.514e-06
## (2 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 7195 1129.9
## 2 7194 1124.6 1 5.3301 34.098 5.469e-09 ***
## 3 7194 1114.1 0 0.0000
## 4 7193 1104.5 1 9.6107 62.590 2.928e-15 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 79416.11
## 2 2a 79384.07
## 3 2b 79316.88
## 4 2c 79256.50
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.08137 0.14102 7.668 1.98e-14 ***
## alpha 0.91493 0.04846 18.879 < 2e-16 ***
## A 159.40185 5.13905 31.018 < 2e-16 ***
## k 35.98938 0.89911 40.028 < 2e-16 ***
## p 0.16407 0.01612 10.178 < 2e-16 ***
## s 2.45792 0.14847 16.555 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.3919 on 7193 degrees of freedom
##
## Number of iterations to convergence: 6
## Achieved convergence tolerance: 9.267e-06
## (2 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 7194 1150.4
## 2 7193 1101.1 1 49.273 321.88 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2c 79256.50
## 2 4 79547.40
## 3 5 79234.25
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.09590 0.14141 7.75 1.05e-14 ***
## alpha 0.91379 0.04844 18.86 < 2e-16 ***
## a 25.29252 2.30219 10.99 < 2e-16 ***
## b 123.73765 4.27083 28.97 < 2e-16 ***
## c 104.02553 3.50472 29.68 < 2e-16 ***
## d 1.33508 0.05096 26.20 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.3913 on 7193 degrees of freedom
##
## Algorithm "port", convergence message: relative convergence (4)
## (2 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 750 92.017
## 2 749 84.644 1 7.3729 65.242 2.646e-15 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 7647.094
## 2 2 7586.205
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.1367 0.2902 0.471 0.638
## alpha 0.9558 0.1105 8.651 < 2e-16 ***
## A 284.8395 36.3007 7.847 1.48e-14 ***
## k 93.3971 16.0340 5.825 8.48e-09 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.3362 on 749 degrees of freedom
##
## Number of iterations to convergence: 5
## Achieved convergence tolerance: 3.457e-07
## (2 observations deleted due to missingness)
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Error in nls(get(paste("f_", Mod.Sel1, "c", sep = "")), data = G_M223, :
## number of iterations exceeded maximum of 50
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 749 84.644
## 2 748 84.115 1 0.52843 4.6991 0.03049 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 7586.205
## 2 2a 7583.489
## 3 2b 6048.497
## 4 2c NA
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau -0.02502 0.27448 -0.091 0.92739
## alpha 0.90237 0.09777 9.230 < 2e-16 ***
## A 214.41085 57.48780 3.730 0.00021 ***
## k 53.39305 23.09154 2.312 0.02111 *
## s 1.22219 0.24916 4.905 1.21e-06 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.3477 on 594 degrees of freedom
##
## Number of iterations to convergence: 9
## Achieved convergence tolerance: 6.551e-06
## (3 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 748 91.700
## 2 747 84.333 1 7.3666 65.251 2.643e-15 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2b 6048.497
## 2 4 7648.498
## 3 5 7587.438
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau -0.02502 0.27448 -0.091 0.92739
## alpha 0.90237 0.09777 9.230 < 2e-16 ***
## A 214.41085 57.48780 3.730 0.00021 ***
## k 53.39305 23.09154 2.312 0.02111 *
## s 1.22219 0.24916 4.905 1.21e-06 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.3477 on 594 degrees of freedom
##
## Number of iterations to convergence: 9
## Achieved convergence tolerance: 6.551e-06
## (3 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 836 173.53
## 2 835 170.37 1 3.1556 15.466 9.099e-05 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 8690.011
## 2 2 8676.613
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.0845 0.6144 1.765 0.0779 .
## alpha 0.6828 0.1673 4.082 4.89e-05 ***
## A 336.2851 72.6004 4.632 4.20e-06 ***
## k 183.2880 41.2191 4.447 9.90e-06 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4517 on 835 degrees of freedom
##
## Number of iterations to convergence: 6
## Achieved convergence tolerance: 8.727e-06
## (1 observation deleted due to missingness)
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 835 170.37
## 2 834 170.22 1 0.1486 0.7280 0.3938
## 3 674 142.26 160 27.9596 0.8279 0.9282
## 4 673 142.26 1 0.0069 0.0327 0.8564
## model AIC
## 1 2 8676.613
## 2 2a 8677.881
## 3 2b 6981.181
## 4 2c 6983.148
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.5623 0.5070 1.109 0.268
## alpha 0.8573 0.1138 7.531 1.63e-13 ***
## A 347.4775 233.5990 1.487 0.137
## k 173.4408 187.9875 0.923 0.357
## s 0.9738 0.1591 6.122 1.57e-09 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4594 on 674 degrees of freedom
##
## Number of iterations to convergence: 6
## Achieved convergence tolerance: 1.693e-06
## (1 observation deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 834 173.48
## 2 833 170.22 1 3.2592 15.949 7.081e-05 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2b 6981.181
## 2 4 8693.783
## 3 5 8679.871
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 0.5623 0.5070 1.109 0.268
## alpha 0.8573 0.1138 7.531 1.63e-13 ***
## A 347.4775 233.5990 1.487 0.137
## k 173.4408 187.9875 0.923 0.357
## s 0.9738 0.1591 6.122 1.57e-09 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4594 on 674 degrees of freedom
##
## Number of iterations to convergence: 6
## Achieved convergence tolerance: 1.693e-06
## (1 observation deleted due to missingness)
## [1] "cannot plot data with prediction"
## [1] "cannot plot observed vs. predicted"
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 1870 1216.8
## 2 1869 1208.6 1 8.2003 12.681 0.0003786 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 23583.09
## 2 2 23572.43
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 5.6743 3.1629 1.794 0.072978 .
## alpha 0.4630 0.1238 3.741 0.000189 ***
## A 242.4795 82.5212 2.938 0.003340 **
## k 231.8338 32.1288 7.216 7.76e-13 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.8041 on 1869 degrees of freedom
##
## Number of iterations to convergence: 6
## Achieved convergence tolerance: 5.895e-06
## (6 observations deleted due to missingness)
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) :
## Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 1869 1208.57
## 2 1868 1192.61 1 15.96 25.0048 6.254e-07 ***
## 3 324 93.71 1544 1098.90 2.4607 < 2.2e-16 ***
## 4 323 93.31 1 0.40 1.3989 0.2378
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 23572.426
## 2 2a 23549.521
## 3 2b 3485.785
## 4 2c 3486.363
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.8945 1.2923 1.466 0.1436
## alpha 0.6580 0.2763 2.381 0.0178 *
## A 65.9248 13.3842 4.926 1.34e-06 ***
## k 30.8188 3.7464 8.226 4.75e-15 ***
## s 42.5340 187.7406 0.227 0.8209
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5378 on 324 degrees of freedom
##
## Number of iterations to convergence: 21
## Achieved convergence tolerance: 8.771e-06
## (1 observation deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 1868 1197.2
## 2 1867 1188.6 1 8.6142 13.531 0.0002413 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2b 3485.785
## 2 4 23556.736
## 3 5 23545.211
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.8945 1.2923 1.466 0.1436
## alpha 0.6580 0.2763 2.381 0.0178 *
## A 65.9248 13.3842 4.926 1.34e-06 ***
## k 30.8188 3.7464 8.226 4.75e-15 ***
## s 42.5340 187.7406 0.227 0.8209
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5378 on 324 degrees of freedom
##
## Number of iterations to convergence: 21
## Achieved convergence tolerance: 8.771e-06
## (1 observation deleted due to missingness)
simple log-normal model: does not fit
log-normal phi model: does not fit
model can fit - but K is negative (only 19 observations) - model excluded
## [1] "cannot plot data with prediction"
## [1] "cannot plot data with prediction"
## [1] "cannot plot data with prediction"
## [1] "cannot plot data with prediction"
## [1] "cannot plot data with prediction"
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 308 124.65
## 2 307 115.58 1 9.0627 24.071 1.509e-06 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 1 3119.284
## 2 2 3097.808
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.6153 1.3279 1.216 0.22476
## alpha 0.8352 0.1541 5.422 1.20e-07 ***
## A 69.6423 17.1618 4.058 6.28e-05 ***
## k 44.9466 15.3373 2.931 0.00364 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.6136 on 307 degrees of freedom
##
## Number of iterations to convergence: 4
## Achieved convergence tolerance: 9.238e-06
## (5 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 307 115.58
## 2 306 115.48 1 0.10112 0.2679 0.6051
## 3 306 114.90 0 0.00000
## 4 305 114.63 1 0.26693 0.7102 0.4000
## model AIC
## 1 2 3097.808
## 2 2a 3099.536
## 3 2b 3097.962
## 4 2c 3099.239
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.6153 1.3279 1.216 0.22476
## alpha 0.8352 0.1541 5.422 1.20e-07 ***
## A 69.6423 17.1618 4.058 6.28e-05 ***
## k 44.9466 15.3373 2.931 0.00364 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.6136 on 307 degrees of freedom
##
## Number of iterations to convergence: 4
## Achieved convergence tolerance: 9.238e-06
## (5 observations deleted due to missingness)
## Analysis of Variance Table
##
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1 306 123.46
## 2 305 115.54 1 7.9259 20.923 6.962e-06 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## model AIC
## 1 2 3097.808
## 2 4 3120.315
## 3 5 3101.680
##
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha *
## B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##
## Parameters:
## Estimate Std. Error t value Pr(>|t|)
## tau 1.6153 1.3279 1.216 0.22476
## alpha 0.8352 0.1541 5.422 1.20e-07 ***
## A 69.6423 17.1618 4.058 6.28e-05 ***
## k 44.9466 15.3373 2.931 0.00364 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.6136 on 307 degrees of freedom
##
## Number of iterations to convergence: 4
## Achieved convergence tolerance: 9.238e-06
## (5 observations deleted due to missingness)
## [1] "cannot plot data with prediction"
| Code | Ecoregion | Sel.Mod.2 | Sel.Mod.3 | Best.Mod |
|---|---|---|---|---|
| 211 | Northeastern Mixed Forest | 2b | 5 | 5 |
| 212 | Laurentian Mixed Forest | 2b | 5 | 5 |
| 221 | Eastern Broadleaf Forest | 2b | 2b | 2b |
| 222 | Midwest Broadleaf Forest | 2a | 5 | 5 |
| 223 | Central Interior Broadleaf Forest | 2c | 2c | 2c |
| 231 | Southeastern Mixed Forest | 2c | 2c | 2c |
| 232 | Outer Coastal Plain Mixed Forest | 2b | 2b | 2b |
| 234 | Lower Mississippi Riverine Forest | 2c | 2c | 2c |
| 242 | Pacific Lowland Mixed Forest | NA | NA | NA |
| 251 | Prairie Parkland (Temperate) | 2b | 5 | 5 |
| 255 | Prairie Parkland (Subtropical) | 2a | 5 | 5 |
| 261 | California Coastal Chaparral Forest and Shrub | NA | NA | NA |
| 262 | California Dry Steppe | NA | NA | NA |
| 263 | California Coastal Steppe - Mixed Forest and Redwood Forest | NA | NA | NA |
| 313 | Colorado Plateau Semi-Desert | NA | NA | NA |
| 315 | Southwest Plateau and Plains Dry Steppe and Shrub | NA | NA | NA |
| 321 | Chihuahuan Semi-Desert | NA | NA | NA |
| 322 | American Semidesert and Desert | NA | NA | NA |
| 331 | Great Plains/Palouse Dry Steppe | NA | NA | NA |
| 332 | Great Plains Steppe | 2b | 2b | 2b |
| 341 | Intermountain Semi-Desert and Desert | NA | NA | NA |
| 342 | Intermountain Semi-Desert | NA | NA | NA |
| 411 | Everglades | NA | NA | NA |
| M211 | Adirondack-New England Mixed forest - Coniferous Forest - Alpine Meadow | 2b | 5 | 5 |
| M221 | Central Appalachian Broadleaf Forest - Coniferous Forest - Meadow | 2c | 5 | 5 |
| M223 | Ozark Broadleaf Forest Meadow | 2b | 2b | 2b |
| M231 | Ouachita Mixed Forest | 2b | 2b | 2b |
| M242 | Cascade Mixed Forest | NA | NA | NA |
| M261 | Sierran Steppe - Mixed Forest - Coniferous Forest - Alpine Meadow | 2b | 2b | 2b |
| M262 | California Coastal Range Coniferous Forest - Open Woodland - Shrub - Meadow | NA | NA | NA |
| M313 | Arizona-New Mexico Mountains Semi-Desert - Open Woodland - Coniferous Forest - Alpine Meadow | NA | NA | NA |
| M331 | Southern Rocky Mountain Steppe - Open Woodland - Coniferous Forest - Alpine Meadow | NA | NA | NA |
| M332 | Middle Rocky Mountain Steppe - Coniferous Forest - Alpine Meadow | NA | NA | NA |
| M333 | Northern Rocky Mountain Steppe - Coniferous Forest - Alpine Meadow | NA | NA | NA |
| M334 | Black Hills Coniferous Forest | 2 | 2 | 2 |
| M341 | Nevada-Utah Mountains Semi-Desert - Coniferous Forest - Alpine Meadow | NA | NA | NA |
| Code | Ecoregion | region | n.obs | n.plots | tau | tau.variance | tau.2.5 | tau.97.5 | alpha | alpha.variance | alpha.2.5 | alpha.97.5 | A | A.2.5 | A.97.5 | k | k.2.5 | k.97.5 | a | a.2.5 | a.97.5 | b | b.se | b.2.5 | b.97.5 | c | c.2.5 | c.97.5 | d | d.2.5 | d.97.5 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 211 | Northeastern Mixed Forest | east | 4838 | 2419 | 1.5826572 | 0.0859229 | 1.0079948 | 2.1573197 | 1.1088849 | 0.0041371 | 0.9827878 | 1.2349819 | 214.23259 | 176.03157 | 252.43361 | 72.38593 | 55.136299 | 89.63556 | 23.88706 | 19.889474 | 27.88464 | 108.49559 | NA | 94.71761 | 122.27358 | 132.60613 | 115.29891 | 149.91335 | 1.194045 | 1.060323 | 1.327767 |
| 212 | Laurentian Mixed Forest | east | 12962 | 6481 | 1.2146588 | 0.0184839 | 0.9481722 | 1.4811455 | 0.8239244 | 0.0013527 | 0.7518336 | 0.8960153 | 112.92651 | 104.53007 | 121.32295 | 41.69107 | 37.451306 | 45.93084 | 13.98018 | 12.400597 | 15.55976 | 71.24841 | NA | 67.14756 | 75.34926 | 116.92499 | 108.58279 | 125.26720 | 1.411748 | 1.327560 | 1.495935 |
| 221 | Eastern Broadleaf Forest | east | 5446 | 2723 | 0.2338300 | 0.0169227 | -0.0211930 | 0.4888529 | 0.8162353 | 0.0009875 | 0.7546312 | 0.8778394 | 293.53650 | 251.49470 | 335.57830 | 64.20469 | 51.624025 | 76.78536 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 222 | Midwest Broadleaf Forest | east | 3552 | 1776 | 1.0026303 | 0.0636318 | 0.5080853 | 1.4971753 | 0.9205052 | 0.0036751 | 0.8016535 | 1.0393569 | 382.71152 | 292.40252 | 473.02052 | 181.08697 | 126.203476 | 235.97047 | 12.37804 | 9.877379 | 14.87869 | 109.05738 | NA | 97.14204 | 120.97273 | 113.94899 | 100.27945 | 127.61854 | 1.184719 | 1.073380 | 1.296058 |
| 223 | Central Interior Broadleaf Forest | east | 6388 | 3194 | 0.8575908 | 0.0148186 | 0.6189556 | 1.0962260 | 0.7120924 | 0.0011431 | 0.6458139 | 0.7783709 | 127.11413 | 118.06312 | 136.16514 | 37.02397 | 34.561601 | 39.48633 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 231 | Southeastern Mixed Forest | east | 7790 | 3895 | 1.5059413 | 0.0337173 | 1.1459915 | 1.8658910 | 0.6511538 | 0.0005098 | 0.6068943 | 0.6954133 | 144.97930 | 132.73182 | 157.22678 | 30.16077 | 27.765806 | 32.55574 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 232 | Outer Coastal Plain Mixed Forest | east | 7940 | 3970 | 0.8355407 | 0.0339820 | 0.4741816 | 1.1968998 | 0.5842890 | 0.0007741 | 0.5297507 | 0.6388272 | 700.20538 | 218.28207 | 1182.12870 | 453.42710 | -103.620712 | 1010.47492 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 234 | Lower Mississippi Riverine Forest | east | 830 | 415 | 0.3953616 | 0.1690247 | -0.4116176 | 1.2023407 | 0.7892546 | 0.0049287 | 0.6514526 | 0.9270566 | 395.47217 | 23.90610 | 767.03825 | 99.96670 | -38.805927 | 238.73932 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 242 | Pacific Lowland Mixed Forest | west | 0 | 0 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 251 | Prairie Parkland (Temperate) | east | 1392 | 696 | 0.2020061 | 0.0639157 | -0.2938224 | 0.6978346 | 0.6558780 | 0.0094616 | 0.4651081 | 0.8466478 | 167.99527 | 131.43342 | 204.55712 | 47.74774 | 34.396282 | 61.09920 | 21.61048 | 16.279939 | 26.94101 | 100.87424 | NA | 86.27605 | 115.47242 | 107.78366 | 88.12748 | 127.43984 | 1.204748 | 1.013088 | 1.396409 |
| 255 | Prairie Parkland (Subtropical) | east | 444 | 222 | -0.0835282 | 0.2320018 | -1.0293926 | 0.8623362 | 0.7987237 | 0.0129602 | 0.5751659 | 1.0222815 | 242.91208 | 126.58518 | 359.23897 | 88.31753 | 28.663971 | 147.97109 | 11.01478 | 6.493934 | 15.53563 | 88.30672 | NA | 67.58342 | 109.03002 | 58.54755 | 48.61408 | 68.48103 | 1.069750 | 0.883204 | 1.256296 |
| 261 | California Coastal Chaparral Forest and Shrub | west | 0 | 0 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 262 | California Dry Steppe | west | 0 | 0 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 263 | California Coastal Steppe - Mixed Forest and Redwood Forest | west | 4 | 2 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 313 | Colorado Plateau Semi-Desert | west | 0 | 0 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 315 | Southwest Plateau and Plains Dry Steppe and Shrub | west | 0 | 0 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 321 | Chihuahuan Semi-Desert | west | 0 | 0 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 322 | American Semidesert and Desert | west | 0 | 0 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 331 | Great Plains/Palouse Dry Steppe | west | 118 | 59 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 332 | Great Plains Steppe | west | 154 | 77 | 2.8327582 | 15.3607851 | -4.9126609 | 10.5781773 | 0.9759876 | 0.0923741 | 0.3753487 | 1.5766264 | 85.62460 | -32.81058 | 204.05978 | 44.97633 | -1.150377 | 91.10303 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 341 | Intermountain Semi-Desert and Desert | west | 4 | 2 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 342 | Intermountain Semi-Desert | west | 2 | 1 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| 411 | Everglades | east | 66 | 33 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| M211 | Adirondack-New England Mixed forest - Coniferous Forest - Alpine Meadow | east | 5108 | 2554 | 1.5624417 | 0.0671905 | 1.0542697 | 2.0706137 | 0.9219623 | 0.0036253 | 0.8039230 | 1.0400017 | 181.37233 | 158.52671 | 204.21795 | 57.00052 | 48.863953 | 65.13709 | 12.85464 | 7.403809 | 18.30547 | 125.02105 | NA | 108.11864 | 141.92346 | 166.06497 | 134.67846 | 197.45149 | 1.546021 | 1.351465 | 1.740576 |
| M221 | Central Appalachian Broadleaf Forest - Coniferous Forest - Meadow | east | 5186 | 2593 | 1.0959008 | 0.0199964 | 0.8186986 | 1.3731030 | 0.9137907 | 0.0023469 | 0.8188253 | 1.0087562 | 159.40185 | 149.32781 | 169.47590 | 35.98938 | 34.226867 | 37.75190 | 25.29252 | 20.779542 | 29.80549 | 123.73765 | NA | 115.36557 | 132.10973 | 104.02553 | 97.15524 | 110.89581 | 1.335080 | 1.235176 | 1.434983 |
| M223 | Ozark Broadleaf Forest Meadow | east | 602 | 301 | -0.0250231 | 0.0753381 | -0.5640881 | 0.5140420 | 0.9023719 | 0.0095586 | 0.7103585 | 1.0943853 | 214.41085 | 101.50679 | 327.31492 | 53.39305 | 8.042058 | 98.74405 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| M231 | Ouachita Mixed Forest | east | 680 | 340 | 0.5622631 | 0.2570724 | -0.4332717 | 1.5577979 | 0.8572806 | 0.0129590 | 0.6337615 | 1.0807997 | 347.47751 | -111.19185 | 806.14687 | 173.44078 | -195.670788 | 542.55234 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| M242 | Cascade Mixed Forest | west | 34 | 17 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| M261 | Sierran Steppe - Mixed Forest - Coniferous Forest - Alpine Meadow | west | 330 | 165 | 1.8945480 | 1.6701111 | -0.6478649 | 4.4369609 | 0.6580466 | 0.0763610 | 0.1144095 | 1.2016837 | 65.92482 | 39.59387 | 92.25577 | 30.81882 | 23.448443 | 38.18920 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| M262 | California Coastal Range Coniferous Forest - Open Woodland - Shrub - Meadow | west | 8 | 4 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| M313 | Arizona-New Mexico Mountains Semi-Desert - Open Woodland - Coniferous Forest - Alpine Meadow | west | 0 | 0 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| M331 | Southern Rocky Mountain Steppe - Open Woodland - Coniferous Forest - Alpine Meadow | west | 0 | 0 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| M332 | Middle Rocky Mountain Steppe - Coniferous Forest - Alpine Meadow | west | 20 | 10 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| M333 | Northern Rocky Mountain Steppe - Coniferous Forest - Alpine Meadow | west | 22 | 11 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| M334 | Black Hills Coniferous Forest | west | 306 | 153 | 1.6153148 | 1.7633679 | -0.9976589 | 4.2282885 | 0.8351977 | 0.0237321 | 0.5320656 | 1.1383298 | 69.64228 | 35.87271 | 103.41185 | 44.94656 | 14.767046 | 75.12608 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
| M341 | Nevada-Utah Mountains Semi-Desert - Coniferous Forest - Alpine Meadow | west | 0 | 0 | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
## OGR data source with driver: ESRI Shapefile
## Source: "C:\Users\hogan.jaaron\Dropbox\FIA_R\Mapping\S_USA.EcoMapProvinces\S_USA.EcoMapProvinces.shp", layer: "S_USA.EcoMapProvinces"
## with 37 features
## It has 17 fields
## Integer64 fields read as strings: PROVINCE_ PROVINCE_I
## Warning: Removed 19 rows containing missing values (`geom_point()`).
## Warning: Removed 19 rows containing missing values (`geom_point()`).